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Abstract 

An explicit quark model, based on a mean field description of non-overlapping 
nucleon bags bound by the self-consistent exchange of a, ui and p mesons, is used 
to investigate the properties of both nuclear and neutron matter. We establish a 
clear understanding of the relationship between this model, which incorporates the 
internal structure of the nucleon, and QHD. Finally, we use the model to study the 
density dependence of the quark condensate in-medium. 
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Recently there has been considerable interest in relativistic calculations of infinite 
nuclear matter as well as dense neutron matter. A relativistic treatment is, of course, 
essential if one aims to deal with the properties of dense matter, including the equation of 
state (EOS). The simplest relativistic model for hadronic matter is the Walecka model [1] 
(often called Quantum Hadrodynamics, i.e., QHD-I [2]), which consists of structureless 
nucleons interacting through the exchange of the a meson and the time component of the 
to meson in the mean-field approximation (MFA). Later Serot and Walecka extended the 
model to incorporate the isovector mesons, tt and p (QHD-II [2, 3]) and used it to discuss 
systems like neutron stars with N ^ Z. QHD has proven to be a very powerful tool for 
the treatment of a wide range of nuclear phenomena [2]. However, it should be noted 
that a strong excitation of nucleon-antinucleon pairs due to the meson fields is somewhat 
unlikely from the point of view of the quark model because the meson fields create only 
one quark-antiquark pair at a time. An extension of such a model to include the structure 
of the nucleon may provide some new insight into the properties of hadronic matter. 

A quark-meson coupling (QMC) mechanism for the saturation of nuclear matter was 
initially proposed by Guichon [4]. In his simple model, nuclear matter consists of non- 
overlapping nucleon bags bound by the self-consistent exchange of scalar (<r) and vector 
(lo) mesons in the MFA. This model has been generalised to include nucleon Fermi motion 
as well as the center-of-mass (cm.) correction to the bag energy [5]. Using this model, it 
has proven possible to investigate both the change of nucleon properties in nuclear matter 
and nuclear structure functions [6, 7]. (An alternative version of this model, based on the 
color dielectric soliton model, has been developed by Naar and Birse [8], who also applied 
it to deep-inelastic scattering.) To treat asymmetric nuclear matter (i.e., N ^ Z) the 
model requires some extension. In this work we show that it is possible to include the 
contribution of the isovector, vector meson, p } in the study of such systems. We shall see 
that it is very interesting to compare this new model with QHD-II and to investigate the 
properties of hadronic matter in terms of quark degrees of freedom. 

Let the mean-field values for the <r, u (the time component) and p (the time component 
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in the third direction of isospin) fields, in uniformly distributed nuclear matter with 
N ^ Z, be cr, uj and 6, respectively. The nucleon is described by the static spherical MIT 
bag [9] in which quarks interact with those mean fields. The Dirac equation for a quark 
field, i/ig, in a bag is then given by 1 



,-0, 



where V a = g q a (T, = g^co and V p = g q p b with the quark-meson coupling constants, g q al g^ 
and g q . The bare quark mass is denoted by m q and T q is the third component of the Pauli 
matrix. (Here we deal with u and d quarks (q=u or d) only.) The normalized, ground 
state for a quark in the nucleon is then given by 

t/nf 3o(xr/R) \ 



v z/3 g <T • fj 1 (xr/R) J V47T 



where 



e q = ttg + R(Vu ±^V P ), for quark 
M~ 2 = 2R 3 j 2 (x)[n q (n q - 1) + Rm;/2]/x 2 



fi g — Rm* 
\\tt q + Rm* ' 



(3) 
(4) 

(5) 



wi 



ith fi„ 



x 2 + (Rm*) 2 and Xg the quark spinor. The effective quark mass, m* is 



defined by 



m* = m q - V a = m q - g q a a. (6) 
The boundary condition at the surface provides the equation for the eigenvalue x: 



JO(X) = f3 q Jl{x). 



(7) 



Using the SU(6) spin-flavor nucleon wave function, the nucleon energy is given by 



E 



N 



P 



E bag + 3V w + ^V p J2r; 



p 



E bag + W u ± -V p , for 



pi 



proton 
. neutron , 



(9) 



1 The sign of the a field value in this paper is opposite to that in Ref.[7] in order to simplify the 
comparison with QHD. 



Here the bag energy is 

E bag = E ^~"° + \*BR\ (10) 

with B the bag constant and z the usual parameter which accounts lor zero-point motion. 

In order to correct lor spurious cm. motion in the bag [10] the mass ol the nucleon 
at rest is taken to be 

Mn = y/E^-ipU), (11) 

where (p 2 cm ) = J2 q {p^} an d (p 2 q ) 1S the expectation value ol the quark momentum squared, 
(x I R) 2 . Inside nuclear matter the effective nucleon mass, Af^, is also given by minimizing 
eq.(ll) with respect to R. To test the sensitivity ol our results to the radius ol the tree 
nucleon, i? , we vary B and z to fix i?o(= 0.6, 0.8, 1.0/ra) and the tree nucleon mass 
(Mjv = 939 MeV). The values ol B and z lor several choices ol radius lor the tree nucleon 



Table 1: B 1 ^ and z lor some bag radii (m q = MeV). 



Ro(fm) 


0.6 


0.8 


1.0 


B^ 4 (MeV) 


188.1 


157.5 


136.3 


zo 


2.030 


1.628 


1.153 



are summarized in table 1. 

The Pauli principle induces Fermi motion ol the nucleon. Since the energy ol a moving 
bag with momentum k is 

1 T . r ( proton 
— Vp, lor 

2 \ neutron, 



e(k) = ^]M% + F + W w ± -V p , lor ( P j ) (12) 



the total energy per nucleon at the nuclear density, is given by 

E tot = — ^( + f kFn )dk^M^ + k 2 + 3V w + ^ ) +-L {m l* 2 -mlu 2 - m y). 

(13) 
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Here p 3 is the difference between the proton and neutron densities, p p — p n} and kp p and 
kp n are the Fermi momenta for protons and neutrons, respectively: 

Pp = T~3 k %> P n = T~^ k Fn and Pb = P P + Pn- (14) 

The to field created by uniformly distributed nucleons is determined by baryon number 
conservation to be 

V w 3gipB 9loPb 



gl mi 



(15) 



(where we define a new coupling constant g w = 3g^) while the a and p mean-fields are 
given by the thermodynamic conditions [2]: 

pfi =° ^ (if) R =»• < i6 » 

\ / R,p B \ / R,p B 



Since the p field value is expressed by 



9p to, (17) 



2m 2 p 



where we take g p = g q } the total energy per nucleon can be rewritten as 



E tot = 2 ,J t Fp + [ kFn )dkJMt 2 + P + ^-p B + ^a 2 + -£—p 3 2 . (18) 
PB^f J J V 2m 2 J 2p B %m 2 p p B V ; 



Therefore, using eq.(16), the value of the a field is given by 

, = /%M^* f^) ■ (19. 



(2tt) 3 7 7 JMtf + P V ^ 



We should note that the expression for the total energy, eq.(18), is identical to that of 
QHD-II. The effect of the internal quark structure enters entirely through the effective 
nucleon mass, eqs.(10) and (11), and the self-consistency condition (SCC) for the a field, 
eq.(19). 

Let us consider the SCC further. The effective quark mass, m*, is given by eq.(6). If 
the nucleon were simply made of three massive constituent quarks the nucleon mass in 
vacuum and that in matter would satisfy 

M N % :!///.. and M* N % 3m* (20) 



and hence, using eq.(6), the effective nucleon mass could be 

M* N = M N - ;t .a. (21) 
where we defined g a = 3g%. Since one finds 

) = ~9a, (22) 

from eq.(21), the SCC gives 

g a a = Zgla (23) 

gl 2 f k Fp f k F r m* 



2(o ,3 v/ +/ ) dk I ^ ( 24 ) 

m l ( 27r ) J J + k 2 

Together with eq.(21), we then find the SCC in the heavy quark mass limit: 

M* = M N -g a a (25) 

m l ( 27r ) J J y/Mff + P 

This is exactly the SCC for QHD-II [2]. 

However, in the more realistic case of light quarks the internal structure of the nucleon 
significantly alters the SCC. In particular, the u and d quarks in the bag model, like current 
quarks, have very small masses and must therefore be treated in a fully relativistic manner. 
Using eqs.(10) and (11), one finds 

M R = ~ 9 ° x Km) [(1 " ^ )S{a) ~ ^ ] (27) 

= -g a x C{a), (28) 
where S(a) is the scalar density of the nucleon bag in matter: 

Hence we find the SCC for the a field in the case where the nucleon has the quark 
structure: 

9l 2 r ,, u f kF P , f k F n . jT M* 



9** = ^7^C(a)( + )dk-==2=. (30) 



Clearly the effect of the internal, quark structure is completely absorbed into the scalar 
density factor, C(a). Thus, we see that it is possible to unify the present QMC model and 
QHD in the MFA: the energy per nucleon and the SCC for the a held are, respectively, 
given by eq.(18) and eq.(19) with 




0.1 0.2 0.3 0.4 0.5 



RV 

Figure 1: Scalar density factors for various bag radii as a function of RV a (m q = 
MeV). The solid, dotted and dashed curves show the results for R = 0.6,0.8 and 1.0/ra, 
respectively. 

three bag radii. We see that it is much smaller than unity and that the dependence on 
the bag radius is not strong. 
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Now we are in a position to present numerical results. First we take m q = MeV, and 
determine the coupling constants, g 2 and g 2 } so as to lit the binding energy ( — 16 MeV) 
and the saturation density (p = 0.17/m~ 3 ) lor equilibrium nuclear matter. Furthermore, 
we choose the p meson coupling constant, g 2 , so as to reproduce the bulk symmetry 
energy, 33.2 MeV [11]. The coupling constants and some calculated properties ol nuclear 

Table 2: Coupling constants and calculated properties ol equilibrium nuclear matter. The 
effective nucleon mass, Af^, the nuclear compressibility, K , and the symmetry energy, a 4 , 
are quoted in MeV. The bottom row is lor QHD-II. We take m a = 550 MeV, m w = 783 
MeV and m = 770 MeV. 







9l/^ 


9 2 P /^ 


M* N 


K 


a 4 


SR 
Ro 


8x 
x 


0.6 


20.2 


1.55 


5.51 


839 


220 


33.2 


-0.03 


-0.07 


0.8 


22.0 


1.14 


5.67 


851 


200 


33.2 


-0.02 


-0.10 


1.0 


22.5 


0.95 


5.69 


856 


190 


33.2 


-0.01 


-0.12 


QHD-II 


7.29 


10.8 


2.93 


522 


540 


33.6 







matter at saturation density are listed in table 2. Because the nucleon is extended these 
coupling constants are related to the coupling constants at the meson poles, g a {im a ), 

Qiuor pi^TTlLuor p) ? by 

g a (im a ) =g a J dfe tk -^ q iP q \^ =m ^ (32) 
gu,, P ( im <o, P ) = 9^,p J dre lk - ? i)\i) q \^ =vm ^ p . (33) 

Using eqs.(32) and (33), the value ol the coupling constant at each meson pole is ob- 
tained: lor example, </^(zm (7 )/47r = 9.35, ^(im a ,)/47r = 5.78 and g 2 (im p ) /Att = 27.4 lor 
Ro = 0.8/ra. It should be pointed out that the strength ol the u field is considerably 
less than that required in QHD because the cm. correction to the bag provides a new 
source ol repulsion [5, 7]. Furthermore, the present model provides values lor the nuclear 
compressibility, K } well within the experimental range [12]: K = 180 ~ 325 MeV. In 
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the last two columns of table 2 we show the relative modifications (with respect to their 
values at zero density) of the bag radius, i?, and the lowest eigenvalue, x, at the saturation 
density. The changes are small (see also Ref. [7]). 




0.5 1 1.5 2 2.5 3 

Pb'Po 

Figure 2: Energy per nucleon for symmetric nuclear matter and for neutron matter. The 
dot-dashed curve is the saturation curve for nuclear matter with R = 0.8/ra. The solid 
curve (for g 2 = 71.2) and the dotted curve (for g 2 = 0) show the results for neutron 
matter with the same bag radius. 

The total energy per nucleon for both symmetric nuclear matter and for neutron matter 
(with Rq = 0.8 fm) are presented in fig. 2. We see that ignoring the p meson coupling yields 
a shallow bound state of neutron matter around ps ~ 0.6/?o 5 but it becomes unbound when 
the p meson contribution is introduced. In fig. 3, we show the equation of state (EOS) 
for neutron matter for both R = 0.6 and 1.0 fm. The p meson contribution is again 
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36 I 1 ■ 1 ' 




12 12.5 13 13.5 14 14.5 15 

log iQ 8 (g/cm 3 ) 

Figure 3: The equation of state for neutron matter. The dashed curve (for R = 0.6/ra) 
and the dotted curve (for R = 1.0 fm) show the results when the rho meson contribution 
is ignored. The solid and dot-dashed curves are the full calculations for R = 0.6 and 
1.0/m, respectively. 

significant in the EOS, while the dependence on the radius of the free nucleon is weak. 
Although the shape of the EOS for neutron matter in the present model is qualitatively 
similar to that in QHDTI, it is much softer than in QHD-II because of the lower nuclear 
compressibility resulting from the internal quark structure of the nucleon. 

Having shown that the QMC model provides an acceptable description of the satu- 
ration properties of nuclear matter we now apply it to study the behaviour of the quark 
condensate in matter. Recent work has linked the quark condensate to a wide range of 
nuclear phenomena [13] and it is an important parameter in any finite density calculation 



10 



using QCD sum-rules [14]. From this point we shall deal only with symmetric nuclear mat- 
ter and introduce a tiny current quark mass lor both u and d quarks: lor example, m q = 10 
MeV. Accordingly, we must readjust the bag parameters and the coupling constants to 
fit the tree nucleon mass and the saturation energy and density. We find that B and z 
become slightly reduced and enhanced, respectively, Irom those values shown in table 1 
and that the coupling constants are 5^/47r = 19.1,20.6,20.8 and g^/^ir = 1.74,1.32,1.13 
lor i? = 0.6,0.8, 1.0/ra, respectively. 

The quark condensate at low density can be related to the nucleon a term, ctjv, which 
is obtained by using the Hellmann-Feynman theorem [15, 16]. Since we find 



dM* N \ 
Jm q ) R 



3C(<r), (34) 



Irom eqs.(6) and (28), the a term lor a nucleon in tree space is 

a N = 3m g C(0), (35) 

where C(a) is the scalar density lactor lor a massive quark. We note that it does not 
change much Irom that in the massless quark case and that the bag radius dependence 
is not strong - unlike fig. 1. II we use C(0) ~ 0.4 and m q = 10 MeV, we obtain ctjv — 12 
MeV, which is much smaller than the empirical value [17]: (T e ^ p ~ 45 MeV. However it 
should be noted that there are many other potential contributions to ctjv, such as the 
meson cloud of the nucleon and its strange quark content [15, 18]. 

Following Cohen et al. [16], the difference between the quark condensates in matter 
and vacuum is given by 

<?<«>- w> = (36) 

1 f d£ dM^ d£ dm a d£ dm^A 
dm a dm q dm^ dm q J ' 

d£ d£ \ 
+ + , (38) 



2 \dM^ dm q dm a dm q dm w dm q 
<j n ( C{a) d£ d£ d£ 

2m q { C(0) dM* N + Xa dm a + Xw d 



where £ is the total energy density, pBE tot} Q(pb) is the in-medium quark condensate, 
(Pb \qq\pb}^ an( l Q(0) that in vacuum. We have used eqs.(34) and (35) and some relations 
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suggested by Cohen et al. [16], = Xa and = Xuj (^) with Xa ~ 



and Xu 



to derive the last equation. Using eqs.(18) and (30) and the Gell-Mann- 



Oakes-Renner relation, 2m q Q(0) = —m^f% } we hnd 



Q{Pb) 
Q(0) 



~ 1 



ml f~ 



m 



9a PO 



[9a (7 



~\2 



m 3 , 



(39) 



where is the pion mass (138 MeV) and /„. the pion decay constant (93 MeV). Note 
that the ratio of the quark condensates, Q(pb) /Q(0), is a simple function of the mean a 
held, <7 CT <7, and ps/po- In fig. 4, the <r held values are plotted at low density. We see that 

300 



250 



200 



1 50 



100 



50 












0.2 



0.4 



0.6 



0.8 



Pb'Po 



Figure 4: Mean held value of the a meson at low densities (with m q = 10 MeV). The 
solid, dotted and dashed curves are for R = 0.6,0.8 and 1.0/ra, respectively. 



the dependence on the radius of the bag is very weak at low density and that the mean 
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field value is well approximated by a linear form for ps < O.bpo: 

g<J a = W a ^ 400 (MeV) {^j . (40) 

Hence, if we assume ctjv = 45 MeV and take g 2 a = 250, g 2 = 18 and C(0) ~ 0.4, we find a 
simple formula for the ratio at low densities 



Q(PB) , n ooi ( Pb\ , ^ f PB 



2N 



w. ^- - :5:51 ^J +O luij- (41) 

which is consistent with the model-independent prediction for the ratio of quark conden- 
sates proposed by Cohen et al. [16]. At ps = 0.5/?o 5 the in-medium quark condensate is 
reduced by roughly 15% from its value in vacuum. 

Furthermore, since the effective nucleon mass can be expanded at low densities as 



(m*) ~ M N (m q ) -Kx^-f 



5 

m*=m Q 



(42) 



L q — " L q 

M N - C(0)( g<T <t), (43) 



we obtain a relation between the quark and nucleon masses 

) , (44) 



Mn ~ m ^ ~_zc(o){- dM * N 

m q — m* I dm* 

where we have used eq.(6). This relation is satisfied within a few percent up to Pb/po 
0.5. If we now solve eq.(43) for Mj^/Mn, we find 



M N ~~ M N K *"~" (45) 
and, using eq.(40) and C(0) ~ 0.4, we find 

#~ 1-0.170 (BS). (46) 
M N \po) y ' 

Comparing this with eq.(41) we see that Q(pb)/Q(0) is essentially (Mfr/M N ) 2 at low 
density. This result is intermediate between the cubic dependence found by Brown and 
Rho [19] and the linear dependence found by Cohen et al. [16]. If on the other hand 
we were to take the limit of heavy quarks considered earlier it is easy to show that one 
recovers the linear dependence. 
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Finally, we would like to add some caveats concerning the present model: (1) The basic 
idea is that the mesons are locally coupled to the quarks. This is certainly common for 
pions in models like the chiral or cloudy bag [20], but may be less justified for vector mesons 
which are not collective states. (2) The model is best suited to low and moderate density 
regions - certainly no higher than ps ~ 3.0/Jo _ because short-range q-q correlations are 
not taken into account. At high densities these would be expected to dominate. (3) The 
pionic cloud of the nucleon [20] should be considered explicitly in any truly quantitative 
study of the properties of the nucleon in-medium. 

In summary, we have presented the quark-meson coupling model based on a mean field 
description of nucleon bags interacting through the exchange of <r, u and p mesons, and 
shown that it can provide an excellent description of the properties of both symmetric and 
asymmetric nuclear matter. The relationship to QHD, and particularly the role of the 
internal structure of the nucleon, was clarified. Finally, the model was shown to provide 
a reasonable description of the quark condensate in-medium. The present model seems 
to be the simplest possible extension of QHD to incorporate explicit quark degrees of 
freedom. 

This work was supported by the Australian Research Council. 
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